Specific heat measurements of the gap structure of the organic superconductors 
K-(ET) 2 Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 
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We present high resolution heat capacity measurements for the organic superconductors 
ft-(ET) 2 Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 in fields up to 14 T. We use the high field data 
to determine the normal state specific heat and hence extract the behavior of the electronic specific 
heat C e i in the superconducting state in zero and finite fields. We find that in both materials for 
T/T c < 0.3, C e i(H = 0) ~ T 2 indicating d -wave superconductivity. Our data are inconsistent with 
s-wave behavior, but may be fitted to a strong coupling d-wave model over the full temperature 
range. 
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The organic superconductors, k— (ET) 2 X, have many 
similarities to the high temperature cuprate supercon- 
ductors (HTSC) Q. In both cases, the electronic 
structure is quasi-two-dimensional and the superconduct- 
ing phase emerges from an antiferromagnetic insulating 
(AFI) state as the phase diagram is transversed. In 
the case of the cuprates, the structure is 'tuned' away 
from the AFI state either by varying the oxygen con- 
tent or making non-isovalent substitutions, whereas in 
k— (ET) 2 X it is achieved either by changing the anion X 
or by applying external pressure. It is natural then to 
speculate that mechanism for superconductivity in these 
two materials may be related even though T c is up to one 
order of magnitude higher in HTSC. 

A first step towards determining if this is indeed the 
case is to determine the symmetry of the superconduct- 
ing energy gap functions in the two families of materials. 
The case of the cuprates has been very well studied and 
the overwhelming consensus is that these materials have 
a gap with predominately d x 2_ y 2 symmetry [ij. In the 
organic materials the situation is more controversial 0]. 

The two most widely studied materials are 
K-(ET) 2 Cu[N(CN) 2 ]Br and k-(ET) 2 Cu(NCS) 2 (here- 
after abbreviated to K-Br and k-NCS) as these have 
the highest T c at ambient pressure (~12K and ~9.5K 
respectively). As yet, no direct phase sensitive de- 
terminations of the gap function have been reported, 
however, there have been numerous experiments which 
have probed its anisotropy. Some early experimental 
determinations of the temperature dependence of the 
penetration depth A(T) reported behavior which was 
consistent with a d-wave gap whereas others were 
consistent with s-wave 0]. The discrepancies stem from 
experiments not having been carried out at sufficiently 
low temperature and/or with high enough precision to 
be conclusive. More recent measurements performed 
down to T/T c < 0.03 showed clearly the existence of low 
energy excitations which were consistent with a <i-wave 
gap [J,@]. Thermal conductivity k(T) data shows a low 
temperature T-linear term [6j and a four fold variation 
with basal plane angle in an applied magnetic field 



0. Both of these results indicate a d-w&ve gap as do 
tunnelling Q and NMR experiments 

Specific heat C has an advantage over many other tech- 
niques in that it is a bulk thermodynamic probe. For 
example, it is largely insensitive to surface contamina- 
tion (or thin layers of damaged material) which may ad- 
versely affect probes such as A(T) (in the Meissner state) 
or tunnelling. A disadvantage is that the electronic com- 
ponent C e i is often only a few percent of the total at 
T c . It is difficult to accurately extract C e / from the 
total which is dominated by phonon contributions. In 
Ref. (loj C e i of K-Br was determined by subtracting an 
estimate of the phonon contribution obtained from a non- 
superconducting quench cooled deuterated version of the 
same compound. These authors found that C e i ~ T 2 as 



However, this approach was 
ll| and Midler et al. [TT 
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expected for a d-wave gap. 
criticized by Elsinger et al. 
instead determined the phonon contribution by applying 
a high magnetic field to destroy the superconductivity. 
These authors claimed that their data for both k-Bi and 
k-NCS were well described by an s-wave gap. It should 
be mentioned however, that in neither of these two re- 
ports were attempts made to explicitly fit their data to 
a d-wave model. 

Here we report high resolution measurements of C e i 
for K-Br and k-NCS which are well described by a strong 
coupling d-wave form of the superconducting gap. The 
low temperature behavior is inconsistent with an s-wave 
gap. 

Samples of both compounds were grown by the usual 
electrochemical method 13j in Argonne, and had masses 
in the range of 80-600 fig. Specific heat measurements 
were conducted in a purpose built calorimeter which uses 
a long relaxation method similar to that described in 
Ref. [3]. Briefly, a Cernox [lj| chip resistor (CX-1030- 
Br) is suspended by silver coated glass fibers in vac- 
uum. The Cernox material acts as both thermometer 
and heater. The sample was attached to the calorimeter 
chip with Apiezon N grease. The addenda (chip, grease, 
and leads) was determined in a separate run immediately 
prior to the main experiment. The thermometer was cal- 
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FIG. 1: (Color online) Left panel: Zero field specific heat data 
for K-Br (abbreviated Br) and «>NCS (abbreviated NCS). Up- 
per right Panel: C/T versus T 2 for for both compounds in 
fields of and 14T. The upper curve for each field is k-NCS. 
Lower right panel: Field dependence of 7 for both compounds. 



ibrated in field, in IT increments up to 14T, by stabi- 
lizing the temperature with a capacitance thermometer. 
The performance of the experiment was extensively stud- 
ied by measurement of high purity samples of Ag (with 
masses in the range 0.3-5mg). In the range 1.3K to 20 K 
the absolute values of the Ag data agreed with standard 
values to within 1%. Small (<4mK) adjustments to the 
calibration points were made to ensure that the measured 
C for all Ag samples were smooth and field independent 
within experimental error. 

Heat capacity data for both materials is shown in 
Fig. [1] In this raw data the ~3% anomaly at T c is barely 
discernable. In order to subtract the large phonon con- 
tribution we have made measurements in magnetic fields 
up to 14T, applied perpendicular to the basal plane. For 
both materials, the maximum field is significantly in ex- 
cess of the upper critical field H C 2, and so at 14T both 
materials are in the normal state. In Fig. [1] we show 
the low temperature portion of the 14T data plotted as 
C/T versus T 2 . Fitting this with a second order polyno- 
mial, we determine the Sommerfeld coefficient 7 as well 
as the coefficients of the leading phonon terms, (3$ and 
/3 5 (C = 7T + f3 3 T 3 + /3 5 T 5 ). The field dependence of 7 
is shown in Fig. Q] and is seen to saturate at 7 = 28 ± 2 
mJ/mol K 2 for > 8T in K-Br and 7 = 35±2 mJ/mol 
K 2 for {j,qH > 3T in k-NCS. In what follows we make the 
assumption that the 14T data is equivalent to that of the 
normal state in zero field (i.e., the only field dependence 
in C is due to the superconductivity). In principle, there 
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FIG. 2: (Color online) Top panel: AC = C(0) - C(14T) 
versus T for K-Br, along with several fits. The thin solid line 
is a fit to the strong coupling d-wave model, whereas the thick 
solid line is the same fit convoluted with a Gaussian. Similarly 
the dashed lines are fits to the strong coupling s-wave model. 
Bottom panel: Enlarged view of the low temperature part of 
the upper panel (the convoluted fits are indistinguishable and 
are omitted for clarity). 



could be magnetic contributions which vary with field. 
However, the insensitivity of C to H at high fields indi- 
cates that these contributions are negligible. Indeed, in 
Ref. [l6l | a sizable magnetic contribution in high field was 
only found for T well below IK. 

In Figs. and 1 we show AC = (7(0) - C(14T) = 
C i — 7T for n-Br and k-NCS respectively. The super- 
conducting transitions are now clearly visible, with the 
midpoint of the transition giving T c —1 2.25 K and 9.56 K 
respectively. Note that the noise in the data is larger 
at higher temperature because the fractional resolution 
of the calorimeter AC/C is roughly constant with T 
whereas the phonon background increases ~ T 3 . 

In many superconductors, the weak coupling form of 
the BCS theory is inadequate to describe in detail the 
physical properties. A full solution to the strong cou- 
pling theory is complicated, and dependent on micro- 
scopic details, but it is found that many properties can be 
explained satisfactorily with the so-called a- model 17 1. 
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FIG. 3: (Color online) The same plots as Fig.[2]but for k-NCS 

Here, the temperature dependence of the energy gap A 
is approximated by the weak-coupling behavior but the 
value at zero temperature is an adjustable parameter. 
This model has been used to describe a wide range of su- 
perconductors including 'exotic' materials such as MgB2 
[H and NbSe 2 0. 

Within the a-model, the entropy S in the supercon- 
ducting state for a two dimensional cylindrical Fermi sur- 
face, is given by 

~^F = — / /m/+[l-/]ln[l-/]&# (1) 

In-Lc T Jo JO 

where the Fermi function / = [cxp^/fcgT) + the 
quasiparticle energy E 2 = e 2 + A 2 (0), j n is the normal 
state 7 and the energy gap A is function of the in-plane 
angle <f>. The specific heat C e i — T§J^. For conven- 
tional isotropic s-wave superconductivity the gap func- 
tion A(cj>, T) = aA^ cs (T), whereas in the simplest case 
for d-wave A(0) = a Ag cg (T) cos 2<fi. In these expres- 
sions Agpg takes the usual s or d-wave weak coupling 
form. 

To allow for the possibility of any part of the sample 
being non-superconducting (and metallic) we allow j n to 
vary in the fit, so the free parameters are a, 7„ and T c . 
As can be seen in the top panels of Figs. [5] and [3J at 



TABLE I: Parameters derived from the s and d wave fits to 
the data in Figs. [2] and [3] The units of 7 are mj/mol K 2 . 714 
is the value of 7„ derived from a fit to the 14T data. The 
maximum gaps at zero temperature Ao = 2.14afc_gT c for the 
d-wave fits. 
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28±2 


1.73 26.6 


1.47 20.0 


rc-NCS 9.56K 


35±2 


1.45 33.3 
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high temperatures the fits to the s and d models are vir- 
tually indistinguishable and both fit the data very well. 
Close to T c the superconducting transition is broadened 
by inhomogeneity and fluctuation effects and the fit is 
considerably improved by convolution with a Gaussian 
(of width a = 0.65 K and 0.43 K for K-Br and k-NCS 
respectively) . 

At lower temperature (bottom panels of Figs. and 
[3]) there is a very significant difference between the two 
models. The <i-wave model fits the data almost exactly 
over the full temperature range whereas the s-wave model 
completely fails at low temperature. The parameters de- 
rived from the fits are given in Table|TJ In both materials, 
for the d-wave fit 7„ is found to be very close to the value 
found from the direct fit to the 14T data 714, whereas the 
s-wave fit is about 30% smaller. If 7„ were fixed at 714 
the s-wave fit would be considerably worse. Clearly, op- 
timizing 7„ for a fit to the low temperature data will 
not significantly improve the s-wave fit. The values of a, 
found show that the d-wave coupling is rather strong. 

The difference between the s-wave and d-wave fits is 
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FIG. 4: (Color online) AC/T versus (T/T c ) for K-Br and 
k-NCS. The solid lines are fits to the d-wave model and the 
dashed and dotted lines are the s-wave fits for each compound 
respectively. 
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perhaps shown more clearly in Fig. [4] where we have 
plotted AC/T versus T/T c . In the low temperature limit, 
the clean <i-wave model predicts C e i ~ T 2 so we expect 
AC/T ~ aT - 7n . The s-wave model predicts AC/T ~ 
a'T~5 exp(— — 7«- The lines on the figure are the 
same fits as in Figs. [2] and [3] Fig. 2] shows that below 
T/T c ~ 0.3 AC/T varies linearly with T, and the full d- 
wave model fits the data over the full temperature range. 
Again, clearly the s-wave model does not fit the at low 
temperature. A linear fit to the AC/T for T/T c < 0.3 
gives a = 2.33 ± 0.03 mJ/mol K 3 and a = 4.21 ± 0.04 
mJ/mol K 3 for K-Br and k-NCS respectively. 

The data in this paper is representative of results taken 
on a large number of different crystals. In total 6 samples 
of K-Br (with T c values in the range 11.5K< T c <12.4K) 
and 3 samples of k-NCS (with T c values in the range 
9.3K< T c <9.6K) were measured and all were found to 
have the same behavior as that reported here. For K-Br 
it is known that fast coolin g th rough the temperature 
region 60-85K depresses T c [201 ] . The sample reported 
here was cooled very slowly (at 0.72 K/hr) through this 
region. Data was also taken for higher cooling rates, 
which we find significantly decreases 7„ and T c but leaves 
the T dependence of C unchanged, except close to T c . 
These results will be reported in detail separately. 

Within the d-w&ve model, by linearizing the angle de- 
pendence of the gap near the nodes it can be shown 
that at low temperature, for a two dimensional cylin- 
drical Fermi surface, ^fjf = 54 ^ 2 3 ^^, 7 " ' where /iA = 
d ^f^ ■ Hence, the coefficient of the T 2 term in the 

^ node 

specific heat determines only the slope of the energy gap 
near the nodes, [iAq. In general, the energy gap may not 
simply vary like cos(2</>) and so /i may differ from 2. We 
have considered this possibility by fitting the data to a 
linearized gap model, where A(<f>) — fJ,Ao4> for \<f>— f | < ^ 
and A(cf>) = A otherwise. The behavior of the data at 
T/T c > 0.3 constrains the values of ^ = 2.0 ± 0.4, al- 
though clearly this depends, to some extent, on our as- 
sumed weak-coupling form of A(T). 

We have also considered the possibility of a mixed or- 
der parameter where the gap does not go to zero at the 
'nodes', for example, d x 2 +y 2 + its. Our data constrains 
e < 0.1, so the minimum gap at the 'nodes' is ^10% of 
the maximum. 

The low temperature (0.3 < T < 2 K) C ei data of 
Nakazawa et aL[l0j also follow a T 2 power law with the 
coefficient of the T 2 term very similar to that reported 
here. However, for T > 2 K their data deviate markedly 
from ours presumably because of the above mentioned 
uncertainties in their phonon subtraction. These uncer- 
tainties are much less important at low temperature. 

Arai et al. Q report that their tunnelling spectra in 
k-NCS can be fitted with a d-wave gap with Ao in the 
range 3.0 — 5.7meV. The lower end of this range compares 
reasonably well to that found here (Aq = 2.6 meV, see 



Table |J). 

Le et al. 5| have reported measurements of the tem- 
perature dependence of the superfluid density by the 
muon spin relaxation (/iSR) technique. We have fit- 
ted their data to the same d-wave model used above 
(again assuming a cylindrical Fermi surface), and find 
a = 1.7 ± 0.2 and a = 1.4 ± 0.2 for re-Br and k-NCS re- 
spectively. This compares favorable to the values found 
here. A quantitative analysis of the penetration depth 
data of Ref. [3] is complicated by uncertainties in the 
assumed value of the penetration depth at zero temper- 
ature, however the measured temperature dependence of 
A is in very good agreement with predictions from the 
current model, once impurity effects are included. 

In conclusion, we have measured the specific heat of an 
extensive set of samples of K-Br and k-NCS and find that 
in all cases the data are well fitted by a strong coupling d- 
wave model. Our data firmly rule out an isotropic s-wave 
gap in these samples. 
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